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Abstract 

We reexamine the role of electron binding effects in the inelastic neutrino- 
atom scattering induced by the neutrino magnetic moment. The differential 
cross section of the process is presented as a sum of the longitudinal and 
transverse components, according to whether the force that the neutrino 
magnetic moment exerts on electrons is parallel or perpendicular to momen- 
tum transfer. The atomic electrons are treated nonrelativistically. On this 
basis, the recent theoretical predictions concerning the magnetic neutrino- 
impact ionization of atoms are critically discussed. Numerical calculations 
are performed for ionization of a hydrogenlike Ge +31 ion by neutrino impact. 

Keywords: neutrino magnetic moment, neutrino-impact ionization 



1. Introduction 

Neutrino electromagnetic properties are among the most intriguing issues 
in neutrino physics (see Ref. [H for a review on this subject). Within the 
Standard Model the value of the neutrino magnetic moment in units of the 
Bohr magneton /i B = e/(2m e ) is 0, 0, S] 
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where m e and m v are the electron and neutrino masses, respectively. There- 
fore, experimental evidence for a value which is much larger than 10~ 19 will 
signalize new physics. 

The nowaday laboratory searches for neutrino magnetic moments are 
dealing with reactor (anti) neutrino sources (MUNU 0, TEXONO Q, and 
GEMMA [3]). The best upper limit for \i v obtained so far in the reactor ex- 
periments is [i v < 3.2 x 10 -11 [t| (see also references in the review article 
which still exceeds by an order of magnitude the most stringent astrophysical 
constraint /i„<3x 10~ 12 jsj]. Since the differential cross section (DCS) for 
the magnetic scattering of a neutrino with energy E u on a free electron (FE) 
at energy transfer T is given by jsl, fiof 



dT m 2 e \T E, 

where a is the fine-structure constant, a way to enhance the sensitivities of 
the experiments is to reduce the low-energy threshold of the detectors with 
respect to the deposited energy T. Note in this connection that the DCS 
due to the weak interaction, d<j( w )/dT, is practically constant in T when 



e v >t m 



The electrons of the detector material are bound to atoms and there- 
fore the FE picture is applicable only when T 3> £&, where Ef, is the electron 
binding energy. If T ~ £&, the electron binding effects must be taken into con- 
sideration. Recently, it was claimed that at T ~ the DCS of the inelastic 
magnetic neutrino scattering on atomic electrons (AE) strongly increases, as 



compared to the FE case (pj), owing to an atomic ionization effect [TlJ. This 
means that the sensitivity of the reactor experiments searching for [i v can 
be significantly enhanced. As demonstrated below, the authors of Ref. [l~D ] 
incorrectly evaluated the DCS, groundlessly employing the approximation of 
real photons in the whole range of the final neutrino angles. In this respect, 
in a more recent work 12] it was argued, by means of a quantum-mechanical 
sum rule, that when E v ^> T and T is comparable to the characteristic 
atomic energies the values of DCS for the AE and FE cases are almost the 
same. While this promising analytical method yields a physically reasonable 
conclusion for the one-electron case, we show in the analysis carried out be- 
low that it requires further elaboration in order to be correctly applied to 
many-electron atomic systems. 

In the present work, we formulate a new theoretical approach which al- 
lows us to consider the magnetic neutrino scattering on atomic electrons in a 
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consistent and physically transparent fashion. This approach can be traced 
down to that developed for the penetration of relativistic charged projectiles 
in matter [13j. According to it, the DCS is given by the sum of two phys- 
ically distinct components-the longitudinal and the transverse, which are 
associated with the corresponding components of the force imposed by the 
neutrino on electrons with respect to the direction of the momentum trans- 
fer q. Such a decomposition also enables us to clearly distinguish between 



the contributions from excited atomic states taken into account in Refs. 11 



and [12 ]. respectively, and to inspect consistently the results of those studies. 
The developed theory is then numerically realized in the case of magnetic 
neutrino-impact ionization (vll) of a hydrogen-like Ge +31 ion. This spe- 
cific case models the .fT-shell ionization of Ge by neutrino impact, which is 
relevant to the reactor experiments 0, 0] ■ 

The article is organized as follows. Sec. [2] delivers general theory for the 
magnetic neutrino scattering on atomic electrons as well as a critical account 



of the theoretical studies published recently [ll|,|l2(. Numerical calculations 



are presented and analyzed in Sec. El The conclusions are drawn in Sec. HJ 
The units h — c — 1 are used throughout unless otherwise stated. 



2. Theory of magnetic neutrino-impact excitation and ionization 
of atoms 

We consider the process where a neutrino with energy E u and momentum 
p„ scatters on an atom at energy-momentum transfer q = (T, q). In what 
follows we focus on the role that plays in this process the interaction of the 
neutrino magnetic moment with atomic electrons, assuming the nucleus to 
be infinitely massive. The latter assumption is reasonable if T S> E^/M, 
where M is the nuclear mass, in which case the recoil of the atom as a whole 
can be safely neglected and the coherent neutrino-atom collisions (without 
changing the atomic state) are of no relevance. The atomic target is supposed 
to be unpolarized and in its ground state |0) with the corresponding energy 
E . We also suppose that T <^ m e and aZ <C 1, where Z is the nuclear 
charge, so that both the initial and the final electronic systems can be treated 
nonrelativistically. The incident and final neutrino states are described by 
the Dirac spinors assuming m u 0. 
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2.1. The differential cross section 

In the low-energy limit the neutrino electromagnetic vertex associated 
with the neutrino magnetic moment has the form 



It should be noted that the \i v related contribution to the neutrino-atom 
scattering does not interfere with that due to weak interaction because in 
contrast to the latter it couples neutrino states with different helicities. Em- 
ploying first-order perturbation theory and using the photon propagator in 
the Coulomb gauge, the transition matrix element for the considered process 
according to Eq. (j2J) is given by 

Mtf = ^(W^Aj((2^-T)(n|f: e ^iO} 



3=1 



+ 2q *'^ wEv-"io>r (3) 



3=1 



where u p \ is the spinor amplitude of the neutrino state with momentum 
p and helicity A, aj the current operator of the jth electron, and \n) the 
final atomic state. Using Eq. (j3J), the double differential cross section can be 
presented as 



dTdQ \dTdQ J i, \ dTdQ J 
where Q = q 2 , T 2 < Q < (2E V - T) 2 , and 

( d 2 a M \ = Tray 2 (2E V - Tf f T 2 



\dTdQj,, m\ AEIQ \} Q 



x 

n 



J2\{n\p(-ci)\0)\ 2 5(T-E n + E f (5) 



d 2 a M \ Tra 2 ^l (2E U - T) 

dTdQ) ± ~ m\ AEIQ 



Q 



{2E V - Tf 

x J]|He ± .j(-q)|0)| 2 5(r -E n + E Q ). (6) 
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In Eqs. (J5]) and (E]), p(— q) and j(— q) are the Fourier transforms of the 
electron density and current density operators, respectively, 



p(-q) = $>*™ (7) 

3=1 



jN-^E^^i + V^)- (8) 

6 3=1 

and the unit vector e_i_ is directed along the p u component which is perpen- 
dicular to q (e_i_ • q = 0). The sums in Eqs. ([5]) and fl6]) run over all atomic 
states \n), with E n being their energies, and, since the ground state |0) is 
unpolarized, do not depend on the direction of q. 

The representation (J4j) is similar to that for the case of charged projec- 
tiles [13| . The longitudinal term ((5]) is associated with atomic excitations 
induced by the force that the neutrino magnetic moment imposes on elec- 
trons in the direction parallel to q. And in the case of nonrelativistic projec- 
tiles, only such excitations are of significance. In the present case, however, 
the projectile is ultrarelativistic. The transverse term ([H]) corresponds to the 
exchange of a virtual photon which is polarized as a real one, that is, per- 
pendicular to q. It can be noted that it resembles a photoabsorption process 
when Q — > T 2 and the virtual-photon four-momentum thus approaches a 
physical value, q 2 — > 0. An important feature is that, due to selections rules, 
Eqs. © and (EJ) involve different, nonintersecting sets of atomic states (see 
Ref. |13j for detail). For this reason they do not interfere. Another marked 
difference between the terms ([5]) and ([6]) is that the first vanishes in the for- 
ward direction, when Q = T 2 , while the second in the backward direction, 
when Q = (2E U — T) 2 . 

In the reactor experiments one typically measures the DCS which derives 
from (jlj) upon integrating over Q: 

da M _ [ {2Ev ' T? d2a M dQ = f <HA + f <HA (9) 



dT J T 2 dTdQ \ dT J ,, \ dT 

Below we separately analyze the longitudinal and transverse components of 
Eq. OU) assuming E u ^> T, which situation is usual for the reactor experi- 
ments. 
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2.2. The longitudinal component of the differential cross section 

When integrating Eq. (J5J) over Q, we neglect the contribution of the term 



oc T 2 /Q that brings only a minor correction to the integral, so that 12 



(^),-^L s(T <' (10) 



where we have introduced the so-called dynamical structure factor |li 

S(T, Q) = |(n|p(-q)|0)| 2 (5(T - E n + E ). (11) 

n 

The latter is related to the density-density or polarization Green's function 
F as follows: 

S(T,Q) = -lmF(T + E ,Q), 

71 

F(T + E ,Q) = (0| p (q) _— i— _ p(-q)|Q) 

= V l(^|p(-q)|0)| 2 (u) 

n 

where H is the atomic Hamiltonian. Using definitions ( II ip and f fT2|) . it is 
straightforward to obtain the dispersion relation 

Fft^lf^^g. (13) 
7T 7o E — E' — iO y ' 

The dynamical structure factor fTTTj) can be expressed in terms of the gener- 



alized oscillator strengths f n [15 



s(t,q) = ^^^(g^r-K + Eo), 

n 

UQ) = ^(E n -E )\(n\p(-cO\0)\ 2 , (14) 



and thus it satisfies the Bethe sum rule 15, 16 



dT 2 7 

S(T,Q)— = —. (15) 
Q m e 
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However, it is not of much help for performing the integration over Q in 
Eq. (TIU|) . For the latter purpose the following momentum-transfer dispersion 
relation was formulated in Ref. [12| (cf. Eq. 

F{E ' Q) -nJ Q>-Q-i dQ > (16) 

provided the energy E is above the ionization threshold. Consider the limit 
Q — > 0. The electron density operator ([7]) at q = is by definition p(0) = Z 
and hence 

F(T + E ,0) = (0|p(0) 1 - p(0)|0) = ^. 

Using it in Eq. (Tl6|) when Q = 0, we arrive at the sum rule (cf. Eq. (1150 ) 

S(T,Q)^- = ^. (17) 



Following the method of Ref. [121 ]. which implies the use of Eq. (fl7|) for 
evaluating the integral in Eq. (fTOl) under assumptions of small T and large 
we get 

da M \ _ 7ia 2 fi 2 u Z 2 _ Z 2<Hl (lg) 



dT ) I, ml T dT 



e 



where daJ^/dT is given by Eq. ([T]) in the limit T/E v — )■ 0. 

Thus, using the dispersion relation ( II 6p . one arrives at Eq. ( lisp where the 
factor of Z 2 occurs, which means that the atomic effects result in a coherent 
enhancement of the DCS as compared to the case of Z free electrons, where 
a typical incoherent-scattering factor of Z is encountered (the same as, for 
instance, in the Compton scattering). This indicates that Eq. f[T6"|) is not 
directly applicable to atoms with more than one electron. At the same time, 
under some modification the method of Ref. 12| can be viable, for instance, 
within the mean-field model of many-electron atomic systems, where the 
occupied one-electron orbitals independently contribute to the vll process. 
In such a case, one employs the sum rule (IT?]) for each electron individually 
(the pro ving of one-electron momentum-transfer sum rules is suggested in 
Ref. [12|). Clearly, this procedure ensures a physically-grounded factor of Z 
in Eq. ( JT8|) after collecting together the one-electron contributions. Note that 
the electron correlation effects beyond the mean-field approximation crucially 
influencing multiple excitation and ionization processes usually slightly affect 
those of single ionization. 



7 



2.3. The transverse component of the differential cross section 

As demonstrated below, it is the case which is relevant to the work 
Integrating Eq. over Q and dropping the insignificant term oc Q/(2E U — 
T) 2 , we get 



dT J m 2 
where the function R is defined as 



1/ R(t,Q)^-, (19) 



R(T, Q) = \(n\e± ■ j(-q)|0)| 2 5(T - £ n + £ ). (20) 

n 

Similarly to the dynamical structure factor ( ITT]) , it can be related to the 
current- current Green's function L (cf. Eq. ( Tl2l ): 

R(T,Q) = -lmL(T + E ,Q), 

L(T + E ,Q) = (0|j ± (q) 1 - - j±(-q)|Q) 

i + i^o — -n _ ^0 

l(«|jx(-q)|0)l 2 (21) 



E 



T - E n + E - iO ' 



where j±(±q) = e_L • j(±q). 

Since j = (pv + vp)/2, with v being the electron velocity operator, we can 
roughly estimate the ratio of the functions ( 120]) and ( fTTl) as ~ t; 2 , where t; a 
1 is a characteristic velocity of atomic electrons. A more accurate estimate 
can be obtained when Q is small on the atomic scale. If Q 2m E Ek one 



can treat Eqs. (TTTj) and (!20]) in the dipole approximation (see also Ref. [121]). 
which yields 

i?(T,Q) T 2 , . 

v = — . (22) 
S(T,Q) Q 1 J 

Note that this ratio is much smaller than unity practically for all Q values 
involved in Eqs. ([TO]) and (JT9"]) . Thus, taking into account the foregoing 
arguments, one might expect the transverse component to play a minor role 



in Eq. ([9]). The authors of Ref. [U|], however, came to the contrary conclusion 
that this component dramatically enhances due to atomic ionization when 
T ~ 6b- The enhancement mechanism proposed in Ref. [llj is based on an 
analogy with the photoionization process. As mentioned above, when Q — > 



8 



T 2 the virtual-photon momentum approaches the physical regime q 2 = 0. In 
this case, we have for the integrand in Eq. (I19p 



R(T,Q) 



Q 



Q ^ T2 An 2 aT 



^ {T) (23) 



where cx 7 (T) is the photoionization cross section at the photon energy T 17 



The limiting form (1231) was used in Ref. ll| in the whole integration interval. 
Such a procedure is obviously incorrect, for the integrand rapidly falls down 
as Q ranges from T 2 up to 4E 2 , especially when Q > r~ 2 , where r a is a char- 



acteristic atomic size (within the Thomas- Fermi model r^ 1 = Z l ^am e 18] ). 
This fact reflects a strong departure from the real-photon regime. For this 
reason we can classify the enhancement of the DCS determined in Ref. [HI 
as spurious. 

3. Numerical estimates of electron binding effects 

Let us inspect numerically the magnetic vll process of a hydrogenlike 
ion with a nuclear charge Z = 32 that mimics the case of a .fT-electron in 
Ge, which is a detector material in the reactor experiments @, 0] • The func- 
tions ( flT|) and ( 120]) are calculated analytically using the well-known expres- 
sions for the generalized oscillator strengths (fT^j) (see, for instance, Ref. fl5[ | ) 



and photoionization matrix elements [17| corresponding to transitions from a 



Is hydrogenlike state. Then, the Q integrations in Eq. are performed nu- 
merically. Fig. [1] shows the DCS (Q, which is normalized to the FE value ([I]), 
above the ionization threshold as a function of T/eb, with the electron bind- 
ing energy given by Eb = a 2 Z 2 m e /2 « 14 keV. As can be seen, the AE results 
for E v ^> 6b are close to the FE ones in magnitude. Such a behavior agrees 



both with the analytical estimates performed in Ref. [12j and with the nu- 



merical results of more rigorous, relativistic treatments for various atomic 



targets 19]. It can be qualitatively explained by noticing the following facts. 
First, in an attractive Coulomb potential there is an infinite set of bound 
states, with the discrete spectrum smoothly transforming into the contin- 
uum at the ionization threshold. Second, the average value of the .fT-electron 
momentum is p e = aZm e and the average change in the electron momentum 
after ejection, Ap e , is such that Ap 2 = 2m e T\ which is analogous to the FE 
case. Thus, in contrast to the claim of Ref. [ll|, at small recoil-electron ener- 
gies one might expect no enhancement of the sensitivity of the experiments 
searching for fi v . 
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Figure 1: The DCS normalized to the FE result versus T/sb for different values of E v . 
The E v = 50eh and E v — lOOe;, curves are practically indistinguishable. 
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Figure 2: The dependence of the ratio between the transverse and longitudinal components 
of the DCS © on T/e b at different values of E v . 
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Fig. [2] presents numerical results for the ratio of the transverse and lon- 
gitudinal terms of Eq. in the same kinematical situations as in Fig. [TJ 
As anticipated, when E u ^> Eb the relative role of the transverse component 
is subsidiary. Though the results for E v = 2eb are not strictly relevant to 
the reactor experiments, where the typical neutrino energies are much higher 
(~ 1 MeV), they are shown for illustrative purposes-in this case the contri- 
bution from the transverse excitations becomes considerable, especially as T 
increases. This finding can be explained by Eq. ( 122]) and different kinematical 
factors 1 - T 2 /Q and 1 - Q/{2E U - T) 2 in Eqs. © and ©, respectively. It 
also can be noted that the kinematical regime where T is close to E v validates 
applicability of the real-photon approximation ( 12 3 p in the whole kinemati- 
cally available region T 2 < Q < (2E U — T) 2 . However, as clearly seen from 
Fig. [H even in a such marginal situation the DCS behaves oppositely to the 
resonant-enhancement scenario of Ref . [Ill ■ 



4. Conclusions 

To summarize, we have performed a theoretical analysis of the magnetic 
neutrino scattering on atomic electrons. For this purpose we have divided 
the DCS into two components corresponding to the longitudinal and trans- 
verse atomic excitations. Incorrectness of the recent theoretical predictions 
assigning the atomic effects to play a significant role in the magnetic neu- 



trino scattering [11] has been revealed. Numerical calculations of the DCS 
for the magnetic vll process of a hydrogenlike Ge +31 ion have been carried 
out at different neutrino energies. The results exhibited suppression of the 
DCS relative to the FE values, which is slight at high impact energy and 
pronounced when the latter is comparable to the electron binding energy. 

The present results based on nonrelativistic calculations qualitatively 
agree with those using a relativistic Dirac-Hartree-Fock treatment of atomic 



electrons [19] . No enhancement of the DCS due to electron binding effects has 



been determined, in contrast to Ref. In this regard, the analytical ap 



proach arguing the insignificance of the atomic effects [12] might have a rich 
potential, provided that it is appropriately modified. Finally, it is unreason- 
able to expect the effects of atomic excitation and/or ionization to introduce 
any enhancement of the sensitivities of the experiments searching for neu- 
trino magnetic moments. Therefore, it will be interesting to explore the role 
of coherent magnetic neutrino scattering on atoms in detectors, which case, 
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however, requires much lower energy thresholds in the deposited energy T 
(~ 100 eV) than presently attainable in the detectors (~ 1 keV). 
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